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Abstract 

Probabilistic interpretation of transition from the dispersive transport regime to the quasi- 
Gaussian one in disordered semiconductors is given in terms of truncated Levy distributions. 
Corresponding transport equations with fractional order derivatives are derived. We discuss 
physical causes leading to truncated waiting time distributions in the process and describe in- 
fluence of truncation on carrier packet form, transient current curves and frequency dependence 
of conductivity. Theoretical results are in a good agreement with experimental facts. 
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1. INTRODUCTION 

Various transport regimes can be observed in disordered 
semiconductors: normal regime characterized by Gaus- 
sian statistics and standard diffusion equation, and differ- 
ent types of anomalous (non-Gaussian) transport. Among 
anomalous transport regimes, the dispersive transport is 
distinguished especially (Scher & Montroll (1975); Zvya- 
gin (1984); Madan & Shaw (1988)). The dispersive 
diffusion packet has a non-Gaussian form, but even so 
it maintains its shape and its spatial extent depends on 
time. In other words, the process reveals the property 
of self-similarity (Scher & Montroll (1975)). Dispersive 
transient current curves are sufficiently differ from the 
normal ones corresponding to Gaussian transport. Current 
decays according to very stretched law, two power law 
sections are picked out: I(t) oc t~ 1+a for t < ix, and 
I(t) oc t~ 1 ~ a for t > tx- The parameter < a < 1 called 
the dispersion parameter depends on material structure 
and transport mechanism. For some mechanisms, temper- 
ature dependence is observed (Zvyagin (1984)). Scher & 
Montroll (1975) have broken through in understanding of 
probability causes of dispersive transport, a is interpreted 
as an exponent of power law asymptotes in a distribution 
of sojourn times T in localized states (waiting times) . Since 
a < 1, the mean value (r) diverges. This leads to non- 
decreasing relative fluctuations of the number of jumps 
between localized states. As a consequence we have non- 
Gaussian form of a diffusion packet and anomalous law of 
packet widening (A oc t a / 2 ). 

The known approaches (such as the Scher-Montroll 
model (1975), the Arkhipov-Rudenko theory (1982, 1983)) 
connect the anomalous-normal transition in the same ma- 
terial with increase of a to 1 because of changing outer 
conditions (for example, temperature). Here, we consider 
another case, the transition occurs due to changes of scale 
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parameters. For small transient times, i. e. small values 
of sample thickness and/or large values of voltages, nor- 
malized transient current curves are practically universal 
and these curves correspond to dispersive transport. But 
in samples with larger thickness, or at smaller voltages 
(at the same temperatures), transient current curves have 
plateau (Bassler (1993); Tyutnev et al. (2005)), that is 
typical for Gaussian transport. Results in this case are of- 
ten described in frames of the quasi-equilibrium transport 
theory. In other words, at increasing transient times, the 
tendency towards the quasi-Gaussian statistics is observed. 
This tendency can not be explained by change of the 
dispersive parameter a, it is a spatiotemporal scale effect. 

Electron transport in polymers is usually modelled as a 
hopping process. Experimental time-of-flight results for 
polymers come to an agreement with the theory, when the 
Gaussian form of an energy distribution of localized states 
is taken (Bassler (1993)). Transition from dispersive type 
of transport to the quasi-Gaussian one can be explained 
by such form of density of localized states. 

The goal of the present paper is to give a probabilistic 
interpretation of transition from the dispersive transport 
regime to the quasi-Gaussian one due to change of the 
scale parameters characterizing the process. The addi- 
tional question is to determine main features of energetic 
and/or spatial distributions of traps leading to the dis- 
cussed phenomenon. Truncated Levy distributions intro- 
duced by Mantegna & Stanley (1994) play an impor- 
tant role in our interpretation. We derive correspondent 
equations with fractional order derivatives and discuss the 
influence of truncation on transient current curves and 
frequency dependence of conductivity. 

2. DISPERSIVE TRANSPORT AND FRACTIONAL 
DIFFERENTIAL APPROACH 

Dispersive (non-Gaussian) transport (Madan & Shaw 
(1988), Zvyagin (1984)) is observed in many disordered 



materials differ in its microscopical structure. A compar- 
ison of available data suggests the presence of universal 
transport properties unrelated to the detailed atomic and 
molecular structure of a substance. The fractional differ- 
ential approach is often applied to description of dispersive 
transport (Barkai (2001); Uchaikin & Sibatov (2005, 
2008); Sibatov & Uchaikin (2009)). 

The Riemann-Liouville type fractional derivatives 
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f(t')dt' 
r(l -a)dt J (t- t') a 



< a < 1 



were for the first time applied to description of dispersive 
transport by Arkhipov et al. (1983b). The authors ex- 
pressed the relationship between concentrations of free and 
localized carriers through the fractional integral. In later 
papers they use a different approximate relation between 
concentrations of localized and free carriers, which they 
called the master dispersive transport equation. This rela- 
tion is believed to hold for any density of localized states 
and permits expressing results through elementary func- 
tions in the case of an exponential density. The Arkhipov- 
Rudenko master equation leads to a diffusion equation 
with a variable diffusion coefficient and mobility (Arkhip- 
iov ct al. (1983a)). 

On the base of the kinetic trapping-emission equations, 
Tiedje (1984) derived a transport equation for free carrier 
concentration. The inverse Laplace transform of this equa- 
tion is nothing but a fractional differential equation (Siba- 
tov & Uchaikin (2007)). 

Barkai (2001) made use of the fractional Fokker-Planck 
equation to account for transient photocurrent relax- 
ation in amorphous semiconductors. He showed the agree- 
ment between selected results of the fractional differen- 
tial approach and results predicted by the Scher-Montroll 
model (1975). 

Power-law decay of photoluminescence in amorphous semi- 
conductors was described in Refs (Seki et al. (2003, 
2006)) in frame of the generalized random walk model with 
recombination by tunnel radiative transitions. Recombina- 
tion was limited by dispersive diffusion of the carriers. In 
the framework of this model, Seki et al. (2006) compiled 
a fractional differential equation for the first passage time 
distribution density. The recombination rate was found 
using the integral Laplace transform of this equation. 

As shown in Refs (Uchaikin (1999); Sibatov & Uchaikin 
(2007)), the main asymptotic terms in the Scher-Montroll 
model satisfy fractional differential equations, the Green 
functions of which are fractionally stable densities. 



3. THE SCHER-MONTROLL MODEL 

Continuous time random walk (CTRW) model, introduced 
by Scher & Montroll (1975), provided the first detailed 
explanation of all the main patterns of current behavior 
in time-of-flight experiments with amorphous semiconduc- 
tors. 

The main assumptions of this model are as follows: 



(1) The transport of charge carriers is a jumplike random 
walk in which the walkers change their positions at 
random instants of time. 

(2) Carrier jumps are independent of one another, and 
time intervals between them (waiting times) are inde- 
pendent, identically distributed random variables T. 

(3) Waiting times are characterized by asymptotically 
power-law distribution: 

P{T > t} otT a , t^oo, 0<a<l. (1) 

Scher & Montroll (1975) simulated charge transfer in 
disordered semiconductors as carrier hopping within the 
model grid of localized states. The grid constitutes a 
regular cubic lattice, each cell of which contains randomly 
distributed sites (localized states). The waiting time till 
the next hopping depends on the distance to the nearest 
neighbor sites. The cell residence time distribution can 
obey the power law owing to site spatial disorder in a cell. 

As is known, the description of normal transport is based 
on the central limit theorem of the probability theory. For 
random quantities distributed according to asymptotically 
power law, divergence of dispersion for a < 2 and diver- 
gence of mathematical expectation for a < 1 make this 
theorem inapplicable, which necessitates the application 
of the generalized limit theorem. 



The generalized limit theorem. Let random quantities 
Xj be independent and identically distributed, and satisfy 
the following conditions 

P(X > x) <~ a + x~ a , x^oo, 
P(X < -x) — a-x~ a , x^oo, 
< a < 2, a + > 0, a_ > and a + + a_ > 0. Then, A n 
and B n > 0, sequences exist such that, as n — > oo, one 
finds 




where S^ a '^ is the stable random variable with exponent 
a and asymmetry parameter /3 = (a + — a_)/(a+ +a_). 
Stable random quantities can be defined through charac- 
teristic functions having the form (form A) (Uchaikin & 
Zolotarev (1999)): 



(k) = cxp{-|fc|"[l-i/3tan(a7r/2) sign(fc)]}, 
5 ( 1 ^(fc)=cxp(-|fc|). 



Certainly, there are an infinite number of sequences 
of normalizing coefficients A n ,B n showing similar 
asymptotic behavior as n — > oo. By way of example, 
they can be defined in the following way (a = (X) and 

c = a + + a_): 



a = 2 A n 

a 6 (1,2) A n 

a = l A n 

a 6 (0, 1) A n 



na, 
na, 

/3cn\nn, 
0, 



B n = Vcnlnn, 

B n = (7rcn/[2r(a)sin(a7r/2)]) 1/a , 
B n = TTcn/2, 

B n = (7rcn/[2r(a)sin(a7r/2)]) 1 / a . 



In the Scher-Montroll model, waiting times T (positive 
random quantities) are distributed according to an asymp- 
totically power law. Therefore, at macroscopic scales, the 
first passage time distribution, conduction current density, 
and concentration of delocalized carriers considered as 



functions of time should have the form of stable density 
distribution. 



The main characteristics of the CTRW model are waiting 
time and jump vector distributions, ip(t) and p(r), re- 
spectively. Spatial distribution density p(r, t) of a particle 
executing random walks and initially located at the origin 
of coordinates is defined in terms of the Fourier-Laplace 
transform (Montroll k Weiss (1965)) 

1 - 



where 



p(k,s) 



p(k,s) 



1 



(2) 



dr e 



ikr 



dt e~ st p(r,t) 



R 

is the Fourier-Laplace transform of normalized particle 
concentration, p(k) is the Fourier transform of path dis- 
tribution density, and ip( s ) is the Laplace transform of 
waiting time distribution density. Substituting into Eq. (2) 
the asymptotic series expansion of the Laplace image of 
waiting time distribution density with the power-law tail 

$(s) ~ l-s a /c a , s<c, 
along with asymptotic expansion of the Fourier image of 
path distribution density: 

p(k) ~ l + ciik-c 2 fc 2 , |k|<l/|ci|, 
and applying the Tauberian theorem, we obtain 

^ „—a r ,a— 1 

p(k, s) = 



„a-l 



-ciik + c 2 k 2 + s a /c a ' 
Rewriting the last expression in the form 

[s a - Cic Q ik + c 2 c a k 2 } p(k, s) = s 
and applying inverse Fourier and Laplace transformations 
yield 

D«p(r, t) + KVp(r, t) - CV 2 p(r, t) - * - S(r), (3) 

1 (1 — a) 

where K and C are vectorial and scalar constants. 
4. TRUNCATED WAITING TIME DISTRIBUTIONS 

It has been emphasized above that the self-similar disper- 
sive transport in disordered semiconductors is character- 
ized by asymptotically power law distributions of sojourn 
times T of carriers (electrons and/or holes) in localized 
states: P{T > t} cx t~ a , t — !> oo, where < a < 1 
is the dispersion parameter. Mean value of such random 
variables diverges. 

It is naturally to suppose that an asymptotically power- 
law distribution of waiting times can be truncated. This 
truncation can be caused by finite values of mobility gap 
at multiple trapping or by limitation of jump lengths 
at hopping. Secondary mechanism acting in parallel to 
the main transport mechanism can be responsible for the 
truncation. We shall consider an influence of truncation 
of power law distributions of waiting times on properties 
of dispersive transport. This influence should become 
apparent in scale effects. 

Mantegna and Stanley (1994) introduced truncated Levy 
flights, a process showing a slow convergence to a Gaus- 
sian. The truncated Levy flight is a Markovian jump 
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Figure 1. Waiting time PDF's for different types of densi- 
ties of localized states in the multiple trapping model: 
exponential (1), Gaussian (2) and rectangular (3). 

process, with the length of jumps showing a power-law 
behavior up to some large scale. At large scales distri- 
bution has cutoffs and thus have finite moments of any 
order. Smoothly (exponentially) truncated Levy flights, 
introduced by Koponen (1995), constructed on Mantegna 
and Stanley's ideas, allows to give a convenient analytic 
representation of results. 

In our model, jump lengths are distributed exponentially 
and waiting times have asymptotically truncated power 
law distributions. We take this distribution in the form 



*(t) = P{T > t} 



7 < c, 



(ct)~ a 
r(l-a) 



exp(-7t), 



(4) 



t > c 



Since the " heavy" power law tail is truncated, all moments 
of random variable are finite, and consequently, a distri- 
bution of a sum of large number of such random variables 
tends to the Gaussian law. However, this convergence is 
very slow and stable Levy distributions play a role of 
intermediate asymptotics. 

In Fig. 1, normalized waiting time distributions numeri- 
cally calculated for three forms of density of states (expo- 
nential, Gaussian and rectangular) in the multiple trap- 
ping model are presented. For all three cases, an average 
energy of states Sq is the same. For exponential DOS 
waiting times are distributes according to the asymptotical 
power law with the exponent a = So/kT. For the case of 
the Gaussian DOS, sojourn times in traps are distributed 
according to a wide law but all moments finite. For the 
case of rectanguler DOS, waiting times have distribution 
that can be approximated by a stretched exponential law. 

Transport is dispersive in all time scales only in the case 
of ideal exponential energy distribution of localized states. 
For two other DOS, transport takes features of the normal 
transport in asymptotics of large times. 



5. TRANSPORT EQUATION FOR THE CASE OF 
TRUNCATED WAITING TIME DISTRIBUTIONS 

Connection between concentrations of delocalized (con- 
duction) n c and trapped nt carriers is expressed by the 
formula 

t 

n t (r,t) = [ T _1 n c (r,*') $(t-t') dt', 



where To is a mean time of delocalized state. Passing on 
to Laplace transforms, taking into account the relation 



e- st *(f) dt~c- a {s + 1 ) t 



s <C c, 



we obtain n t (r, s) ~ c~ Q T ~ 1 (s + 7) a ~ 1 n c (r, s). The inverse 
Laplace transformation leads to the equation with deriva- 
tive of fractional order: 



n c (r,t)~c a T e-^ D 



\ 1 — a ft 



t > c 



(5) 



Substituting the last relation into the continuity equation 
<9n(r, t) 



dt 



div ^fiE n c (r, t) — D Vn c (r, i) 







and taking into account the fact that most of carriers are 
trapped, i. e. n(r,t) w nt(r,t), we arrive at 
<9n(r, t) 
fit + 

+div [e' 1 \D 1 - a e' lt (K n(r, t) - C Vn(r, t))] = 0, (6) 

where K = c Q ro//E, C = c a ToD are coefhcients, and 
H, D, E are a mobility, a diffusion coefficient and a field 
intensity, respectively. If a — > 1, Eq. (6) becomes the 
standard Fokker-Planck equation describing the normal 
transport. When 7 = 0, the equation coincides with the 
dispersive transport equation. 

6. SCALE EFFECT OF TRANSITION FROM THE 
DISPERSIVE REGIME TO THE GAUSSIAN ONE 

Fig. 2 represents numerically calculated spatial distribu- 
tions for one-sided random walks. With different waiting 
time distributions, one can see the transition from frac- 
tionally stable statistics to Gaussian statistics for waiting 
time distribution having truncated power law tail. The up- 
per graph represents results for exponentially distributed 
times between jumps. In the middle graph, waiting times 
are not truncated, the exponent of the power law tail is 
equal to 0.5. In the lower graph, waiting time distribution 
has the sharp truncation at time r tr = 1500 (arbitrary 
units). The exponential distribution of waiting times has 
the same mathematical expectation as truncated power 
law in the lower graph. In the case of exponentially dis- 
tributed waiting times, we see the fast convergence to the 
Gaussian distribution, in the second case, the distribution 
becomes fractionally stable law at some time t and main- 
tains the form for all following times. In case of truncated 
power law tails, the crossover between fractionally stable 
and Gaussian statistics is observed. 

Conduction current density at pulsed injection for the 
case of exponentially truncated power-law waiting time 
distributions is expressed as 

j(x,t) = eTVexp [jy-J -7^J X 
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Figure 2. The transition from the fractionally stable statis- 
tics to the Gaussian one in the case of a one-sided 
random walk with waiting times distributed according 
to law having truncated power law tail. Solid lines 
represent Gaussian PDFs, dotted lines correspond to 
fractionally stable densities. 
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Transient current density is found by substituting this 
expression into the formula 



1 



I(t) = - / j(x,t)dx. 



If a = 1/2, the expression for transient current takes the 
form: 



I{f), arb. units 



Re[cj(co)/G(0)] 




Figure 3. Transient photocurrent in the case of truncated 
power law distribution of waiting times for different 
values of L/Z-ratio. 
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Figure 4. Theoretical frequency dependencies of conduc- 
tivity for different a values. 

This formula is important in that "a knowledge of the 
fluctuations of the equilibrium ensemble in the absence 
of the electric field permits a calculation of the linear 
response of the system (mobility)" (Scher & Lax (1973a)). 
Scher & Lax showed that relation (10) can be written out 
as: 



+V7 



erf(VT^) - erf ( y/jt 



2^ 



(8) 



D(u) 



Fig. 3 illustrates transformation of transient current curves 
with increasing L/Z-ratio. When the time of flight is 
much smaller than the truncation time 7 _1 , the trans- 
port remains dispersive and does not pass to Gaussian 
asymptotics. If ttr is compared with r y~ 1 the shape of 
the transient current curves undergoes modification, and 
they become inconsistent with the curves for normal and 
dispersive transport. For i t r 7 _1 , transport in the long- 
time asymptotic regime becomes quasi-Gaussian. 

7. FREQUENCY DEPENDENCE OF 
CONDUCTIVITY 

The frequency dependence of the real component of con- 
ductivity in disordered semiconductors is usually fairly 
well described by the power law: 

Re cr(cj) = Aj 7 , (9) 

where the exponent 7 normally takes on values from 0.7 
to 1 (Zvyagin (1984)). The dependence of such type is 
characteristic of a very broad class of materials. 

Conductivity is related to mobility by the expression 

<j{lo) = er][x(ui). 
Here, 77 is the concentration of effective carriers. The 
Nyquist formula (generalized Einstein relation) linking 
mobility with the diffusion coefficient at nonzero frequen- 
cies has the form 

= (e/kT)D(u>), 
where the noise spectrum according the Wiener-Khintchin 
theorem is expressed through the Fourier transform of the 
velocity autocorrelation function 



Re D(w) = / cos(ut)(v(t)v(0))dt. 



(10) 



1 s 



dt e 



<[r(i)-r(0)] 2 ). (11) 



The latter relation is possible to rewrite in the form 



D(oS) = —uj 2 I dx x 2 \n(x, s)] € 



(12) 



where nix, s) is the Laplace image of n(a;, t) with respect 
to t. 

Equation (6) for one-dimensional case without field as- 
sumes the form: 

dn(x,t) = _^ d 1 -* lt d 2 njx,t) 

dt at 1 -* dx 2 ' 

The Laplace transform of this equation yields 

8 nix, s) =C(s + 7) 1 "" + Six). 

Substituting the solution of this equation, 



~, \ s -l/2( s+7 )(a-l)/2 

n{x, s) = — exp 



tc V Vc 

into relation (12) gives 

D(u) = 2Cij + iuj) 1 - a . 

Hence follows 

Re ct(cj) = (e 2 r;/fcT) Re D{u) = 
= 2C(e 2 r,/kT) ( 7 2 + aj 2 ) (1_Q)/2 cos((l - a) arctan(w/ 7 )). (13) 

For frequencies uj 3> 7, it is easily shown that 

Re ct(cj) = 2Cie 2 rj/kT) w 1 "" sin(vra/2). 



Figure 4 represents frequency-dependent conductivity 
curves calculated from Eq. (13). Formula (13) predicts 
the power-law dependence of conductivity on lo at high 
frequencies in the dispersive transport case. The exponent 
may acquire values from to 1; in normal transport (a — > 



1), conductivity is totally frequency independent. In trans- 
port driven by the multiple trapping mechanism, expo- 
nent a grows linearly with temperature. Consequently, the 
exponent s = 1 — a in the frequency dependence of con- 
ductivity in the case of alternating current must linearly 
decrease with increasing temperature. Such temperature 
behavior has been reported for a variety of semiconductors 
(see, for instance, Ghosh et al. (2006)). 

8. CONCLUSION 

Here, transition from the dispersive transport regime to 
the quasi-Gaussian one in disordered semiconductors is 
interpreted in terms of truncated Levy distributions of 
waiting times. So, polymer with Gaussian density of local- 
ized states is not exclusive representative of materials that 
can show such behavior. The phenomenon is more general 
and it is based on statistical rules such as generalized limit 
theorem. Analytical results are supported by numerical 
simulations. 
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